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Recent experimental demonstrations of quantum coherence of the charge and flux states of Josephson 
junctions show that the quantum Josephson dynamics can be used to develop scalable quantum logic 
circuits. In this work, I review the basic concepts of Josephson tunneling and Josephson-j unction 
qubits, and discuss two problems of this tunneling motivated by quantum computing applications. 
One is the theory of photon-assisted resonant flux tunneling in SQUID systems used to demonstrate 
quantum coherence of flux. Another is the problem of quantum measurement of charge with the 
SET electrometer. It is shown that the SET electrometer at the Coulomb blockade threshold is the 
quantum-limited detector with energy sensitivity reaching h/y/S in the resonant-tunneling regime. 



I. INTRODUCTION 

Remarkable ability of Bose condensates to create a sin- 
gle quantum state occupied with a macroscopically large 
number of particles leads to an unusual situation when a 
macroscopic object behaves as a pure quantum state of a 
simple quantum mechanical system. This behavior is in 
stark contrast to the usual situation when macroscopic 
objects behave purely classically and do not show any 
collective quantum effects despite the underlying quan- 
tum dynamics of their microscopic elements. The macro- 
scopic quantum states of Bose condensates can exhibit 
non-trivial quantum dynamics associated with the fact 
that the number of particles in the condensate and the 
phase ofrits wavefunction are non-commuting quantum 
variablefjilla. 

Macroscopic quantum phenomena arising from this dy- 
namics are best studied in Josephson tunnel junctions, 
where the role of bosons is played by Cooper pairs of 
electrons, and both the number of pairs and the phase 
of the wavefunction give well-defined and controlled con- 
tributions to the junction energy. The "number" part 
of the energy comes from the Coulomb charging of the 
junction capacitance C by Cooper pairs and has charac- 
teristic magnitude Ec = (2e)^/2C, while the pair tun- 
neling across the junction creates the "phase" term in 
the energy of magnitude Ej. The nature of the quan- 
tum dynamics of a Josephson junction depends on the 
ratio of these two characteristic energies. In large junc- 
tions with Ej ^ Ec the phase behaves almost classically 
with quantum effects being "corrections" to semiclassical 
dynamics. One of the first quantum effects studied in 
this regime is the macroscopic quantiup tunneling of the 
phase - see reviews in the collectioncl. In the opposite 
case of small junctions with Ej <C Ec, the junction dy- 
namics can be interpreted as semiclassical dynamics of 
charge on the capacitance C. Characteristic feature of 
the charge dynamics is correlated tunneling of individ- 
ual Cooper pairs one Iw .one between superconducting 
electrodes of the systemQ B. 

Macroscopic character of the quantum effects in 



Josephson junction dynamics makes them attractive as 
the basis for development of scalable quantum logic cir- 
cuits for quantum computation. The regime of coher- 
ent quantum dynamics in Josephson tunneling can be 
reached in relatively large, micrometer-size junctions that 
can be designed and fabricated with large degree of preci- 
sion. The relevant junction parameters, including the en- 
ergies Ej and Ec, reflect the collective properties of the 
superconductors forming the junction and are not sensi- 
tive to their precise microscopic structure on the atomic 
scale. This makes the Josephson junctions much less 
susceptible to the microscopic disorder unavoidable in 
the present-day fabrication technology, than, e.g., semi- 
conductor quantum dots in the regime of coherent elec- 
tron transport, also considered for quantum computing 
applicationa]. 

The advantage of scalability of Josepshon junction and, 
in general, solid-state qubits and logic circuits comes at 
a price of stronger dephasing and energy relaxation in 
these systems than in other non-solid-state qubits. Nev- 
ertheless, recent experimental demonstrations of quan- 
tum nnherent Josephson dynamics both in chargeu and 
phaseBM regime show that the problem of dephasing and 
energy relaxation in Josephson junctions is in principle 
solvable, and stimulate further studies of superconduct- 
ing qubits. The aim of this work is to review the basic 
concepts of quantum dynamics of mesoscopic Josephson 
junctions and qubits, and to discuss two specific ques- 
tions motivated by quantum computing applications of 
Josephson junctions: the theory of coherent effects in the 
photon-assisted resonant flux tunneling, and the problem 
of quantum measurement with SET electrometer. 



II. JOSEPHSON TUNNELING IN MESOSCOPIC 
JOSEPHSON JUNCTIONS 



An isolated mesoscopic Josephson tunnel junction is a 
system of two small superconducting islands separated 
by an insulator layer that is sufficiently thin to allow 
electrons to tunnel between the islands. The junction 
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is "mesoscopic" if its capacitance C is small enough for 
the charging energy Ec = [2ef /2C of a single Cooper 
pair to be sufficiently large. Quantitative description of 
this system is most transparent when the superconduct- 
ing energy gap of the islands is much larger than other 
characteristic junction energies, including charging en- 
ergy Ec and temperature T . In this case, quasiparticles 
can not be excited in the superconductors and the junc- 
tion properties can be understood simply from the stan- 
dard notion that the superconductors are Bose-Einstein 
condensates of Cooper pairs, with Cooper pairs occupy- 
ing a single quantum state in each island. If ni and 
712 are the numbers of Cooper pairs in the two supercon- 
ducting islands, the imbalance n of the number of Cooper 
pairs between the two superconductors, n = (ni — n2)/2, 
gives the charge 2en on the junction capacitance C . The 
junction Hamiltonian depends on n through the charging 
energy (2en)^/2C— 2enVe, where Ve is the potential dif- 
ference induced between the two superconducting islands 
by external electric field. 

Another part of the Hamiltonian comes from the tun- 
neling of Cooper pairs between the islands. In accordance 
with the standard practice, we denote the amplitude of 
this tunneling as —Ej/2 and take into account that the 
transfer of a Cooper pair between the islands changes n 
by ±1. Then, combining the electrostatic and the tunnel- 
ing terms, we get the total Hamiltonian H of the junction 
in the n-representation: 

H = Ecin - qf - ^{\n){n + 1| + \n + l){n\) . (1) 

Here q = CVe/2e is the charge (in units of 2e) induced 
on the junction capacitance by the source of potential 
Ve- The energy Ej of Cooper pair tunneling is usually 
referred to as Josephson coupling energy. It can be ex- 
pressed in terms of the normal-state tunnel resistance of 
the junction. 

Behavior of a single Josephson junction described by 
the Hamiltonian (|l|) is non-trivial due to important dif- 
ference between the two charges n and q that define its 
charging energy. The charge n originates from the total 
number of Cooper pairs in a superconducting island and 
can take on only integer values, as a charge of any isolated 
conductor. By contrast, the charge q is the polarization 
charge induced on the capacitance C by external elec- 
tric field and can be varied continuously. Because of this 
difference, jfiigenstates of the Hamiltonian (|l]) form en- 
ergy bandstj periodic in q with the period one. When q 
changes by one, a Cooper pair is transferred through the 
junction increasing n also by one, and fully compensat- 
ing for the increase of q, whereas smaller changes of q can 
not be compensated by any transfer of integer number of 
Cooper pairs. 

An example of the periodic dependence of junction 
properties on q is given in Fig. ^ which shows one pe- 
riod of g-dependence of the ground state energy Eq of the 
junction obtained by direct numerical diagonalization of 



the Hamiltonian (P for several values of the Ej / Ec ra- 
tio. Periodic oscillations of the ground-state energy, and 
similar oscillations of energies Em of the excited states, 
result in oscillations of all other junction characteristics 
with q. For instance, differentiation of the junction free 
energy F = —Tln"^^^ e^^-^l^ , with repect to q gives the 
average number of Cooper pairs (n) charging the junc- 
tion: 



(n) 



1 dF 
2E^'~d^ 



(2) 



and the average voltage across the junction V{q) 
2e{{n)~q)/C. 
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FIG. 1. The ground state energy (a) and the average num- 
ber of Cooper pairs on the junction capacitance (b) for an 
isolated mesoscopic Josephson junction, as functions of the 
external charge q injected into the junction. The curves in 
(b) are for the same values of the ratio Ec/ Ej as in (a). 

Average number of Cooper pairs (||) on the junction 
capacitance as a function of q is shown in Fig. |^b for zero 
temperature and different Ej/Ec ratios. In the limit 
of small coupling energies £'j, the number of pairs (n) 
increases in a step-wise manner remaining constant for 
almost all of q values exept for half-integer g's, when it 
rapidly increases by one. Qualitatively, at Ec ^ Ej, the 
supercurrent flow between the junction electrodes is "dis- 
cretized" into the transfer of individual Cooper pairs by 
strong Coulomb repulsion, and variations of the charge 
q injected into the junction leads to the possibility of 
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controlling their tunneling. 

Step-like behavior of (n) at Ej ^ Ec reflects quan- 
tization of the number of Cooper pairs in the supercon- 
ducting electrodes of the junction. When Ej increases, 
the junction electrodes become strongly coupled and the 
quantization of n is gradually washed out. In this case, 
the ground state energy of the junction (Fig. [l^) becomes 
virtually independent of q, and the average number (n) of 
Cooper pairs on the junction follows closely the induced 
charge q (Fig. 0b), (n) ~ q. Qualitatively, this suppres- 
sion of n quantization with increasing Ej is the result 
of large quantum fluctuations of n induced by tunneling 
terms in the Hamiltonian (^. 

At Ej > Ec, the junction properties can be described 
more directly by introducing the phase operator ip canon- 
ically conjugated to the number n of Coopet pairs. From 
quantum mechanics of a single boson modeEl, it is known 
that the operator of the phase ipj of the condensate wave- 
function in one of the junction electrodes, j — 1, 2, can 
be introduced through the operators s"^ = \n ± l){n\ 
(omitting the index j to simplify the notations) raising 
and lowering the number of Cooper pairs in the con- 
densate. On all of the states |n) besides the state |0) 
with no particles, the operators satisfy the relations 
s+s^ = s^s+ — 1 and — {s^Y which suggest that 
they can be expressed as 

s± = e±'^ . (3) 

In fact, since the state |0) is annihilated by s^, s^|0) = 0, 
the commutator of s+ and does not vanish: 

[s-,s+] = \0)(0\. (4) 

Non-vanishing commutator contradicts the possibility of 
expressing the operators through one phase operator, 
and leads to the well-known problem with a generic def- 
inition of the phase of a boson mode, discussed iu the 
literature on quantum optics - see, e.g., the reviewli3 and 
references therein. 

However, in the case of a Bose condensate with macro- 
scopically large number of particles, all relevant states 
do not have the component |0), the commutator (^) 
effectively vanishes, and eqs. (^) define the operators 
giiip^ Then, the commutator between n and the rais- 
ing/lowering operators can be written as 

e^"^{n±l) ^ne^"^ . (5) 

For a superconducting island, the number of Cooper pairs 
in the condensate can be estimated as the ratio A/ 6 of 
the superconducting energy gap A over the spacing S 
of the single-particle energy levels. Condition A 3> 5 
necessary for defining the phase operator caiucides with 
the criterion of the BCS superconductivitylljo. 

With the raising and lowering operators e^'*^' for the 
two condensates in the junction electrodes, the opera- 
tor of transfer of a Cooper pair through the junction 
can be written as 6=^*"^ = e^^'^'^^^'^^l They have the 



same commutators (|5|) with the number of Cooper pairs 
n on the junction capacitance C. The operators e^*''' 
allow one to define any periodic function of ip, and the 
commutation relations (^ imply that in this periodic (p- 
representation, the number of Cooper piars can be writ- 
ten as n = —id /dip. (The aperiodic operator of the phase 
Lp itself can not be defined for an isolated junction, since 
n is integer and all junction properties should be peri- 
odic in ip with the period 2tt). The Hamiltonian ([^) is 
equivalent then to that of a plain rotor in an external 
field: 

H = Ec{-^-qf -Ejcos^. (6) 

I Oip 

This representation of the Hamiltonian (|l]) allows one to 
analyze simply the behavior of the ground state energy 
eo and average charge (n) (Fig. 0) for Ej > Ec in terms 
of tunneling through the coupling potential —Ej cos (p. 

Descrpition in terms of the phase dynamics becomes 
particular relevant in the situation of an "rf SQUID" 
(see Fig. ^ below), when the junction is included in ex- 
ternal superconducting loop. Due to quantization of the 
magnetic fiux $ in the loop (see, e.g.,t3) the phase tp is 
linked to the flux, (p = 27r<I>/<I>o, where $o = T^fi/e is the 
"flux quantum" . The states of this system that differ in 
phase by 27r are clearly distinguishable, and the charge 
n = —id/dip is no loger integer. The transition from pe- 
riodic to aperiodic phase ip in the Josephson junctions 
has been debated extensively in the eatb^-discussions of 
charging effects in Josephson junctionsOETllfl. In phys- 
ical terms, in the rf SQUID, both electrodes of the junc- 
tion are parts of the same conductor. In this situation, 
the charge n is the polarization charge and can vary con- 
tinuously. One consequence of this is that any external 
charge q can be screened by charge transfer through the 
loop, so that q does not have any effect on the junc- 
tion dynamics and can be omitted from the Hamilto- 
nian. Then, the Hamiltonian of the system, including 
the magnetic energy of the loop inductance L threaded 
by an external flux $e is: 

H^Ecn^-Ejco.^-,^^^^^!:^^-^. (7) 

In such an "aperiodic" system, the phase operator ip is 
well-defined and has the standard commutation relation 
with the charge nn 

[ip,n]^i. (8) 

Both the controlled quantum dynamics of individual 
Cooper pairs governed by the Hamiltonian (|l|) and quan- 
tum dynamics of flux (R) can be used to build qubits 
for quantum computation. The main obstacle to this 
is dephasing and relaxation mechanisms present in the 
Josephson junction circuits that couple charge or flux 
motion to other degrees of freedom not included in eqs. 
(|l|) or (^. Such a coupling makes dynamics irreversible 
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and suppresses quantum correlations of the charge or flux 
states. Although suppression of the relaxation strength 
to acceptably low level is a challenging problem, re- 
versible quantum coherent dynamics has been demon- 
strated both in the charge and flux regimeS|. In the 
charge case, this has been done implicitlyEjEll, through 
observation of splitting of the two lowest energy eigen- 
states of the Hamiltonian (||) at q ~ 1/2, and explic- 
itly, by direct measurement of the time-dependent oscil- 
lations of a Cooper paiia. Quantum coherence of the flux 
states remained peiusive for many years. Very recently, 
two experimentsEHl3 reported first observations of the co- 
herent superposition of the macroscopically distinct flux 
state, with superposition of states evidenced indirectly 
by the coherent splitting of the energy eigenstates of the 
Hamiltonian similar to (|^). The splitting is caused by 
tunneling through the maximum of the SQUID potential 
[/(<!>) for the flux dynamics. In the next section, I dis- 
cuss briefly the theory of flux tunneling in one of these 
experiments. Mare detailed presentation of this theory 
can be found inE3. 



III. QUANTUM COHERENCE OF FLUX STATES 

With the appropriate choice of parameters of the 
SQUID Hamiltonian (0), $e =i $o/2 and Ej ~ $g/27r2L, 
the SQUID potential U{^) has two minima separated 
by a potential barrier. Quantum coherence of flux state 
should manifest itself as coherent oscillations of flux be- 
tween the two potential wells due to tunneling through 
the barrier. The. evidence for such coherent oscillations 
was obtained inO by studying photon-assisted resonant 
tunneling of flux between the two well. lEhe energy di- 
agram of this process is shown in Fig. ||.tj An rf signal 
of frequency uj resonantly couples the ground state |0) 
in the left well of the potential to an excited state |1) 
with energy E in this well. The condition of the reso- 
nant excitation is that the detuning v = E — lu is small, 
V <^ E. If the amplitude a of the rf excitation is also rel- 
atively small, a <^ E, photon-induced coupling between 
the states |0) and |1) can be described in the rotating- 
wave approximation, in which the terms in the coupling 
Hamiltonian that oscillate rapidly in time are neglected. 
The excited state |1) in the left well is coupled resonantly 
with amplitude — A/2 to a state |2) in the right well that 
is shifted in energy by e with respect to |1). The to- 
tal Hamiltonian for the flux dynamics in the basis of the 
three states |0), |1), and |2) is: 

/ a/2 \ 
H ^ \ a/2 V . (9) 

\ -A/2 v-e j 

The rate of resonant flux tunneling betweea-jthe. 
wells of the SQUID potential in absence of rf signalEfrEEl 
is known to be qualitatively independent of the strength 
of the energy relaxation in the flux dynamics. The shape 



of the resonant tunneling peaks is the same both in the 
regimes of coherent and incoherent tunneling. In contrast 
to this, the rate of the photon-assisted tunneling is sensi- 
tive to the relaxation strength, with the resonant peaks 
acquiring additional structure in the regime of coherent 
tunneling, when the energy relaxation is weak. 
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FIG. 2. Schematic diagram of the "photon-assisted" 
macroscopic resonant tunneling of flux stimulated by an rf 
perturbation of strength a/2. 

The main source of the energy relaxation for quantum 
dynamics of Josephson junctions is electromagnetic en- 
vironment of the junction. Under the assumption that 
electromagnetic modes of the environment are in equilib- 
rium at temperature T, the interaction between the flux 
$ and the heat bath of these modes can be written as 

= -//$, (10) 

where is the fluctuating current created in the SQUID 
loop by the environment, with the correlation function 
given by the fluctuation-dissipation theorem: 

Here the brackets (...) denote averaging over the equilib- 
rium density matrix of the environment, and G(ui) is the 
dissipative part of the environment conductance. 

The coupling ( [lO| ) leads to the two types of relaxation 
processes for dynamics of flux tunneling in the double- 
well potential. One is the "intrawell" relaxation inducing 
transitions from the states |1) and |2) to the states with 
lower energies within the same well. In these transitions, 
the environment is absorbing large energy quanta on the 
order of energy separation between the states in the wells. 
Another is "interwell" relaxation which arises from the 
fluctuating energy shifts of one well relative to another 
caused by the coupling to environment (p^. Interwell 
transitions transfer to the environment much smaller en- 
ergy on the order of energies in the Hamiltonian (^) . 

The dissipative dynamics of flux with the Hamiltonian 
(^) can be described with the time-evolution equations 



4 



for the density matrix p of the system. These equations 
can be obtained in the standard way from the environ- 
ment characteriS|tif s ( p^ ) , (^l|) in the case of weak relax- 
ation - see, e.g.j£3. Under the natural assumption that 
the average flux in the states |0) and |1) in the left well of 
the SQUID potential is the same, the part of the coupling 
( p^ ) that corresponds to the interwell relaxation is: 



i 2P22 



(16) 



V = 




(12) 



where (5$ is half of the difference between the average 
flux in the left and right wells. While eq. (|l^) is written 
in the flux basis |0) , |1) , |2), weak relaxation is conve- 
niently described in the basis of the eigenstates \n) of 
the Hamiltonian (|^). In this basis, the contribution of 
the weak interwell relaxation ( |l2|) to the evolution of the 
density matrix p is given by the standard expressions: 



Pnn — ^ ^ (S/mn Pmra ^nniPnn) ; 



(13) 



Transition and dephasing rates in these equations are: 

— ~ I _ _ JJ \2 

]nm — ^ i Inm ^ '-'mm) j 

where J7„m are the matrix elements of the operator U 
( |l^ ) in the eigenstates basis, e„ are the eigenenergies, and 
dimensionless parameter g = 2G{6^)'^ /Ti characterizes 
the strength of the interwell relaxation. 

The intrawell relaxation has simple form directly in the 
flux basis. Equations for the off-diagonal elements of p 
in this basis are: 

Poi = (ii" - ri/2)/9oi + ia{poQ - Pii)/2 - zApo2/2 , 

P12 = -{is + (Ti + r2)/2)pi2 + iA{p22 - Pii)/2 - iapo2/2 , 

P02 = («(!^ - e) - r2/2)po2 - tapi2/2 - iApoi/2 . (14) 

The rates ri_2 in these equations describe, respectively, 
the relaxation out of the states |1) and |2) into the lower- 
energy states within the left and right wells. Diagonal 
part of the evolution equations is: 



The evolution equations for p that include both the 
interwell and intrawell relaxation can be solved numeri- 
cally. To do this, one diagonalizes the Hamiltonian (||), 
calculates the interwell relaxation terms (^sj) in the eigen- 
states basis and transforms them into the flux basis where 
the the intrawell relaxation has the simple form (|l4|), 
([Tsl). Finally, the stationary value of the density matrix 
p gives the flux tunneling rate ([l^). 

Figures |3| shows an example of the results of such a 
calculation obtained for vanishing temperature T and 
Fi = r2 = F. The tunneling rate is plotted as a func- 
tion of the detuning ly for e = and several values of the 
relative strength of the intrawell relaxation F. We see 
that the main resonant tunneling peaks is split in the two 
peaks due to coherent oscillations of flux between the two 
wells. Appearance of such splitting can be easily under- 
stood, since sufficiently weak rf signal excites the system 
not into the state |1) localized in the left well, but into 
the two hybridized states formed out of the states |1) and 
1 2) in the two wells. The two different energies of the two 
hybridized states lead to the two peaks in the tunneling 
rate. In the time-domain, hybridization of the states |1) 
and 1 2) corresponds to the coherent oscillations of flux 
between these states. This means that the splitting of 
the resonant tunneling peak is the direct manifestation 
of the quantum coherent oscillations of flux between the 
two wells of the SQUID potential. 



Poo = -almpoi + Fipii + V2P22 , 
pii = almpoi + AImpi2 - Vipn , 
P22 = -AImpi2 - F2P22 • 



(15) 



Equation for poo here was written under the assumption 
that after the relaxation in the right well the system is 
returned instantly to the initial state |0). This modifica- 
tion of the evolution equations allows us to calculate the 
rate r^^ of the flux tunneling between the two wells frouji 
the stationary solution p^^"^ of the evolution equationsEa: 



a=0.3A 
e/A=0.0 
g=0.03 




FIG. 3. The rate r~ of the photon assisted flux tunneling 
as a function of the detuning 1/ in the case of symmetric cou- 
pling between the tunneling flux states, e = 0, in presence of 
both the interwell and intrawell relaxation. DiflFerent curves 
correspond to different magnitude of the intrawell relaxation 
rate F relative to the interwell relaxation. 

At £ = 0, the hybridized states are distributed sym- 
metrically between the two wells, and asymmetry of the 
two peaks in Fig. ^ is caused by the interwell relaxation. 
The positive-i^ side of the double-peak structure corre- 
sponds to excitation of the system into the lower-energy 
eigenstate and is unaffected by the interwell relaxation 
at zero temperature, since there is no energy in this 
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regime to create additional tunneling path. In contrast, 
the negative-z/ side of the double-peak structure corre- 
sponds to excitation of the system into the eigenstate 
with larger energy, and the interwell relaxation increases 
the rate of tunneling out of this state. Because of this, 
the negative-z/ peak in Fig. || is larger that the tunnel- 
ing peak at positive ly, and the tunneling rate at < 
decreases much more slowly away from the peak than at 
positive ly. 

In absence of the interwell relaxation, the double-peak 
structure is symmetric for e = as illustrated in Fig. ^. 
This figure also shows how the splitting of the tunneling 
peak is suppressed with increasing relaxation strength. 
The curves are obtained by finding the stationary solu- 
tion of eqs. (|lj) and ( [l5| ) which are valid for arbitrary 
strength of the intrawell relaxation, Fi 2 A, a, and ne- 
glecting the weak interwell relaxation ([l^ ) in this regime. 
They show that for F > A the two peaks of the double- 
peak structure in the flux tunneling rate are suppressed 
and merge into one broad resonant tunneling peak. 
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FIG. 4. Evolution of the double-peak structure in the rate 
of the photon-assisted resonant flux tunneling as a func- 
tion of detuning f with increasing relaxation rate F. 

It should be noted that for sufficiently large bias e, the 
two peaks of the resonant flux tunneling exist even for 
strong relaxation rate F, as can be seen analytically in the 
limit Fi_2 ^ A, a. Expressing the off-diagonal element of 
p in terms of the diagonal ones in the stationary regime 
from eqs. (14) and using their stationary values in eqs. 

( p^ ) we find the stationary probability p^22 ^^'^ 
tunneling rate (iq). For strong relaxation. 
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(4t.2+r2)(4(i.-£)2+F2) 



(17) 



At large bias energy e, the tunneling rate ( p^ ) as a func- 
tion of detuning v contains two separate peaks that cor- 
respond to excitation of the flux into the energy states 
1 1,2) localized in the two wells. The two peaks merge, 
however, with decreasing bias and form one broad peak 
at e ~ shown in Fig. ^ This behavior should be con- 



trasted with the coherent-tunneling regime with small F, 
when the two peaks lemain separated even at e ~ 0. 

In the experimentu, the flux tunneling rate was 
studied as a function of the external flux <I>e through the 
SQUID loop. Variations of $2 change both the energy 
shift £ between the two wells of the SQUID potential 
U{^) and the separation E between the levels within the 
same well. For a fixed frequency u of external rf ra- 
diation, detuning v = E — lj changes together with E. 
Although there are some minor quantitative differences 
between the dependence of the tunneling rate r^^ on $e 
(which controls both v and e), and dependence of 
on the detuning v plotted in Figs. ^ and ^, both depen- 
dences show the splitting of the resonant tunneling peaks 
due to coherent superposition of the flux states in theiwo 
wells of the potential. This means that observational of 
the splitting of the resonant peaks of the photon-assisted 
flux tunneling under the resonant bias conditions e ~ 
demonstrates quantum coherent flux dynamics and quan- 
tum superposition of macroscopically distinct flux states. 
This coherent flux dynamics can be used for development 
of the flux qubits. 



IV. CHARGE AND FLUX QUBITS 

Josephson junction qubits and quantum logic circuits 
can operate in the charge and flux regimes. To create 
a qubit, one needs to reduce the junction dynamics to 
two states and to be able to control accurately the re- 
sulting two-state dynamics. The most basic arrangement 
that satisfy these requipements in the charge space is the 
single-Cooper-pair boxE3"E2l configuration, in which the 
external charge q is injected into the junction by a gate 
voltage through an additional capacitance (Fig. ^). 



U(n) 




FIG. 5. Equivalent circuit (a) and the energy diagram (b) 
of the charge qubit. External gate voltage V coupled through 
a small capacitance Co induces charge close to half a Cooper 
pair charge 2e on the Josephson junction with coupling energy 
Ej and capacitance C. 

The charge qubit is the Cooper-pair box with small 
Josephson coupling energy Ej <C Ec and induced charge 
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q ~ 1/2.E11e3 In this case, the two charge states, n = 
and n = 1, are nearly-degenerate and are separated from 
all other states by the large energy gaps on the order of 
Ec (Fig- ||b). The two states are coupled by the tunnel- 
ing of one Cooper pair, and the junction Hamiltonian (Q) 
is effectively reduced to the two-state Hamiltonian: 

H = Ec{q - i)(|0)(0| - |1)(1|) - ^(|0)(1| + |1)(0|) . 

(18) 

The energy difference between the charge states is con- 
trolled through q by the gate voltage, while their tunnel 
coupling can be controlled by the external magnetic field 
that modulates the Josephson couphng energy Ej. Op- 
eration of a Cooper-pair box as a qubit has been demon- 
strated experimentally intl. 




FIG. 6. Equivalent circuit (a) and the energy diagram (b) 
of the flux qubit. External flux $e induces the flux through 
the inductance L of the rf SQUID that is close to half a flux 
quantum. 

The basic scheme of the flux qubitEHsl (Fig. ||) is sim- 
ilar to the charge qubit. When external flux through an 
rf SQUID (0) is equal to a half of the magnetic flux quan- 
tum $0, and Ej ~ <^1/2tt^L, the SQUID potential U{<^) 
has two symmetric minima at $ = $o/2 ± A$ with the 
two degenerate ground states in them. These two states 
are separated from the excited states by large energy 
gaps on the order of quantum of oscillation frequency ujp 
around the minima, hujp ~ (EjEcY^"^ , and are coupled 
by tunneling through the barrier separating the minima. 
The tunneling amplitude A/2 is small, A/2 ^ hujp. In 
this regime, the SQUID dynamics reduces to the two- 
states dynamics at temperatures T <^ hujp. The en- 
ergy difference e between the two flux states is controlled 
by the deviations of the external flux $e from $o/2, 
£ ~ ($e — $o/2)A$/i, while the tunneling amplitude can 
be controlled by additioaal flux modulating the Joseph- 
son coupling energy Ejc3. This means that if the rate of 
the environment-induced decoherence of the flux states 
is suppressed significantly below the tunneling frequency 
A, as in experiments£Hl3, the rf SQUID can be used as 
the flux qubit. 

After individual qubits, the next step in the design of 



the quantum computer is the organization of the qubit 
dynamics necessary to perform quantum logic operations. 
The main general requirement to this dynamics is that it 
should maintain the level of precision sufficient for quan- 
tum computation. For instance, the amplitudes of the 
unwanted transitions to the higher-energy states of the 
system (leading to leai;age of the probability out of the 
two-state qubit spacd^ and destroying correct operation 
of the qubit) should be low. The best and probably 
the only type of dynamics that satisfies this requirement 
in the solid-state systems is the adiabatic dynamicsEil. 
When all qubit parameters are changed adiabatically, the 
amplitudes of all unwanted transitions decrease exponen- 
tially with the ratio of the characteristic qubit energies 
{Ej for the charge and A for the flux qubits) to the 
frequency of qubit operation. For comparison, the am- 
plitude of the unwanted transitions ip4^he regime of Rabi 
oscillatipBS employed in the ion-trapE3 or NMR qubits - 
see, e.g.LO, decreases only linearly with the operation fre- 
quency/energy ratio. This means that despite the term 
"adiabatic" , for a given level of accuracy and given qubit 
energies, adiabatic manipulations of the qubit states can 
be performed considerably faster than Rabi transitions. 
Larger operation frequency of adiabatic quantum logic 
gates is particularly important for solid-state gates which 
should have significantly shorter decoherence time that 
the NMR or ion-trap gates. 

One more difficult problem of the organization of the 
quantum logic circuits out of individual qubits is the ne- 
cessity to switch on and off the interaction between the 
qubits in the controlled manner. Parameters of the di- 
rect Coulomb interaction between the charge qubits and 
magnetic interaction between the fiux qubits typically de- 
pend on the system geometry and can not be controlled 
by external voltage or current signals. Two possible so- 
lutions of this problem has been suggested so far. One is 
the spatiaLseparation of the qubit states in the arrays of 
junctionsEll, another is making the interaction between 
the qubits dependent on the "inteipal" qubit parame- 
ter and controlling this parameterEj. The advantage of 
the first approach is that it gives the possibility of con- 
trolling the interaction with sufficiently high accuracy, 
as needed for quantum computation. However, the diffi- 
culty of fabricating and controllingJijing arrays of junc- 
tions make simpler approaches liketJ attractive for the 
near-term experiments. 

V. SINGLE-ELECTRON TRANSISTOR AS 
QUANTUM DETECTOR 

Essential part of the operation of a quantum circuit is 
the measurement process. In the case of the quantum 
computing circuits, the measurement is necessary for the 
read-out of the final state, for error-correction during op- 
eration, and, possibly, for preparation of the initial state. 
The standard detector for measurements of charge in sys- 
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tems of mesoscopic tunnel itmctions is the single-electron 
tunneling (SET) transistoiOcHl. |0|V|ejj the last ten years, 
it has been extensively develop edEJea and-pused in many 
experiments as electrometer - see, e.g.,L2rc3c3. In this 
Section, we study the detector properties of SET tran- 
sistor in a quantum measurement. Quantum measure- 
ments with the standard flux detector, a dc SQUID, are 
discussed inEll. j— . 

SET transistorey is a small conductor, i.e., semicon- 
ductor quantum dot or small metallic island, placed be- 
tween two external electrodes (Fig. |^). The conductor 
forms two tunnel junctions with the electrodes and acts 
as the central electrode of the transistor. Due to Coulomb 
charging of this electrode by tunneling electrons, the cur- 
rent I through the junctions can be controlled by the 
electrostatic potential (j) of the central electrode. Sensi- 
tivity of the current / to variations of (j) makes it possible 
to detect motion of charges in the vicinity of the tran- 
sistor, and is the basis for its operation as electrometer. 
When the bias voltage V across the transistor is smaller 
than the threshold Vt of Coulomb blockade, the tunnel- 
ing is suppressed by the charging energy that is required 
to transfer an electron in or out of the central electrode. 
In this regime, the current / flows through the transis- 
tor only by a two-step "macroscopic quantum tunneling" 
(MQT) process via a virtual intermediate-state below the 
Coulomb energy barrier (for review, seeE3), and is much 
smaller than the current above the Coulomb blockade 
threshold Vt- 




FIG. 7. Diagram of the SET transistor as electrometer. 
The tunnel current / (measurement output) indicated by ar- 
rows is controlled by electrostatic potential 4> of the central 
electrode of the transistor (input signal) created by a mea- 
sured system, and is coupled to the transistor through the 
charge Q of its central electrode. 

Rapid current rise aX V ~ Vt makes the Coulomb 
blockade threshold the optimuip-poerating point of the 
SET transistor as electromctcrBo. Motivated by this 
fact, we analyze below characteristics of the SET elec- 
trometer as a quantum detector close to the Coulomb 
blockade threshold. Transport properties of the SET 
transistor in the vicinity of the threshold are determined 
by the crossover from MQT to the regular single-electron 
tunneling. The dominant feature of this crossover is 
quantum broadening of the transistor charge states by 



tunnelingL3e3~Ea. Such a broadening determines the 
shape of the current-voltage characteristics I{V) at F ~ 
Vt , and can be described in the simplest way in the regime 
when the temperature and "bias energy" eV are smaller 
than the energy spacing of the single-particle states of 
the central electrode. The tunnel current / flows then 
only through one single-particle state and the problem 
becomes equivalent to the standard problem of resonant 
tunneling through an individual single-particle level. Al- 
though this is not the regime of operation of a typical Al- 
based SET electrometer, it is routinely realized in trapsr 
port through semiconductor quantuip-|dots - see, e.g.,E3, 
and is reached even in Al transistoralS. 

The resonant tunneling through one single-particle 
level can be described quantitatively with the standard 
Hamiltonian: 

H = ea^a + ^ (£pc]pCjp + Tjpc]pa + T'^a'^Cjp) . (19) 

3P 

The operators Cjp represent electrons in the state p with 
energy £p of the jth electrode of the transistor, j = 1, 2, 
while operator a represents electron occupying the reso- 
nant level. The energy e of the resonant level depends on 
the electrostatic potential of the central electrode of the 
transistor, 5e = —eScj). Potential </> controls the tunnel 
current / through the transistor and serves as the input 
signal of the measurement. We assume that the mea- 
sured system is coupled weakly to the transistor, so that 
the typical variations of and e are small and transistor 
operates as a linear detector. In the resonant-tunneling 
regime this implies that Se <^ eV, 7, where 7 is the tun- 
neling width of the resonant level. We also assume that 
the signal frequency is much smaller than eV and 7 and 
treat the transistors response as instantaneous. 

The characteristics of a linear detector relevant for 
quantum measurement are the output noise, back-action 
noise of the detector on the measured system, and the 
linear response coeflicient of transformation of the input 
into the output signal - see,e.g.,tJ. For the SET tran- 
sistor, the measurement output is the current / through 
the junctions. The output noise is then the current noise 
with spectral density Si{uj), while the response coeffi- 
cient A is equal to 5{I)/d(j), where 5{I) is the change of 
the average current / due to variation 6(1) of the poten- 
tial (j). When the potential is produced by a dynamic 
system, e.g., a qubit, the transistor is coupled to this 
system by coupling energy ipQ, where Q = ea'a is the 
operator of the charge on the resonant level. Since the 
charge Q fluctuates due to the current flow through the 
resonant level, the coupling creates a fluctuating force 
acting on the measured system that leads to dephasing 
of the quantum states of this system. Such a "backac- 
tion dephasing" is the fundamental aspect of a quantum 
measurement that is responsible for localization of the 
measured system in the eigenstates of the measured ob- 
servable. In the case of SET transistor, the measured 
observable is the potential (f>, and intensity of the back- 
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action dephasing is determined by the spectral density of 
the charge noise Sq{uj). 

In the assumed regime of low signal frequency, only the 
zero- frequency values of the spectral densities of both the 
output and the input noise, as well as the spectral density 
Siq{uj) of the correlator between I and Q are relevant for 
the measurement problem. The zero-frequency value of 
S'/q(w) can be calculated as 



SiQ^—Jdt[{I{t)Q)-{I){Q)], (20) 

with the zero- frequency limits of Si{lu) and Sq{uj) given 
by similar expressions. 

Because of the simple quadratic form of the Hamilto- 
nian ( |l9| ) all these characteristics of the transistor can 
be found exactly. Partially, this was done in the the- 
ory of charge sensitivity ottlie SET electrometer in the 
resonant-tunneling regimeEd'Ell. Calculation of the charge 
sensitivity for measurements of static charges requires 
only the knowledge of the current noise S*/, whereas under 
a generic measurement conditions one needs to know also 
the backaction noise Sq and the correlation Siq. ftlow 
all spectral densities are calculated following Ref.Ej. In 
contrast to other treatments of noise properties of the res- 
onant tunneling structuresHEj which focus on the sup- 
pression of shot noise, discussion of the optimum detector 
properties of such structures requires the consideration 
of the noise properties at temperatures and bias voltages 
smaller than the width 7 of the resonant level. 

The Heisenberg equations of motion for operators cjp 
and a obtained from the Haniiltonian (|l9[) are: 



icjp — EpCjp Tjpa , ia — ea 



(21) 



Integrating the first of eqs. ( ^ ) we get: 



Cjp{t) ^ Cjp{t) + iTjp dt'e-*^''(*-*')a(t') , (22) 



where Cjp is the "initial value" of Cjp that describes the 
equilibrium occupation of electron states in the transistor 
electrodes: 



(23) 



Here fj{£p) are the Fermi distribution functions with the 
chemical potentials iij of the two electrodes, fJ-i — fJ-2 — 
eV > 0. 

The second of eqs. (^ij) combined with eq. ( p2|) can be 
transformed as 

/•OO 

a = ^lea + i ^ T*^,p - ^ \T,p\^ / dt' e^'^'' a{t - t') . 

■JO 



(24) 

Under the assumption that the electrodes have continu- 
ous spectrum of energy states, and that the densities of 



states Pj{sp) together with the matrix elements of tun- 
neling are constant in a wide energy range around ~ e, 
the last term in this equation reduces to —ja{t), with 



7 = 7i + 72 , 1] 



Pj{e). 



Here |Tjp is the tunneling probability \Tjp\'^ averaged 
over all states with the same energy Cp, and 7^ is half 
of the tunneling rate between the jth electrode and the 
resonant level. After this transformation, eq. (|3) gives: 



jp 



(25) 



Equation 



allows us to find the backaction noise 



Sq^^ I dt{a''{t)a){a{t)a'< 



■E 

1,3 



111] 



dv 



[{e - i/)2 +72 



,U{y){l-f.^{y)). (26) 



In the limit of low temperatures, T <^ eV,j, relevant 
for transistor operation as a quantum detector, eq. ( [2^ ) 
reduces to the following form: 



Q 



e 7172 



dv 



1^)2 -I- 72]2 



(27) 



When either temperature T or bias eV are larger than 
the tunneling width 7, integral in eq. (|2^) can be taken 
treating /'s as constants, and gives expression for the 
noise Sq that can also be obtained from solution of a 
simple kinetic equation for the occupation probabilities 
of the resonant level. The same conclusion applies to eq. 
( p7| ) in the case of large voltages, eV :$> when the 
noise Sq in the occupation of the resonant level is due 
the shot noise of the tunneling current. In this case, eq. 
gives: 



Sq 



e^7i72 
2777^ 



(28) 



In order to find the average current (/) and the noise 
Si at low frequencies, the tunneling current I through the 
transistor can be calculated in either of the junctions. In 
the first junction. 



*e^(ripc{ 



(29) 



Using explicit expression ( |25| ) for a in eq. ( |22| ) to express 
Cjp in terms of Cjp, we write the operator ( p9| ) of the 
current / as 



ie{A'^a- a'^ A) , 



(30) 



where 



9 



3P 



z- Sp + i^ 



(31) 



Averaging eq. (|30|) over the equilibrium state of the elec- 
trodes as in the eq. (p3|), we get the standard result for 
the do current through a resonant level, and find the re- 
sponse coefiicient A = S{I)/5(j): 



A 



4e^7i72 



dv 



r(/iH-/2H). (32) 



At large bias, when e — fi2^ J and = — £ 7, the 
response coefficient decreases, e.g., for low temperatures, 



A = 



2e^7i72 
ttE^ 



The noise Si can also be found from eq. (|3 
to eqs. ( pO| ) and (^), we get: 



(33) 
Similarly 



Si^^ dt[{A\t)A){a{t)a^) + {A{t)A^){a^ it)a)^ 
zvr 



e{I) 27172 
27r 7^ 



(36) 



where the average current through the level is (/) = 
267172/7. In structures that are not highly asymmet- 
ric, 71 ~ 72, eq. ( p6| ) describes suppression of the cur- 
rent noise in comparison with the regular shot noise 
e(/)/27r. The previous discussion shows that this sup- 
pression has two correct interpretations that appear to 
be quite different. One explanation relates the suppres- 
sion to the quantum Fermi correlations in the tunneling 
flow of electrons through the structure described by the 
transparency-dependent suppression factor (1 — D). In 
another interpretation, the suppression is brought about 
by the classical correlations in the process of electron 
transfer through the structure, which break this process 
into two distinct steps, tunneling onto the resonant level 
and out of it. 

Finally, following the same steps that lead to eq. ( |3^ ) 
for the current noise, we calculate the current-charge cor- 
relator SiQ. From eqs. ( pO| ) and (|30|), and the definition 
of Q we have: 



{A\t)a){a{t)A'') - {A(t)A^) (a\t)a)] . 

Evaluating averages in this expression with the help of 
eqs. (|l]) and (p5|), we get after some algebra: 



e 7172 



dv 



[(e- j/)2 +72 



SiQ^^ dt[{A^it)a){a{t)a^) - {Ait)a^){aHt)a)] . 
Ztt J 

Making use of eqs. ( p5| ) and ( ^ ) to calculate the averages 
in this expression we obtain: 



SiQ ~ 



6^7172 



dv 



[(e - z^)2 -I- 72 



r[7i/i(^)(l-/i(^))^ 



[7i72(/i(^^)(l-/i(;^)) + /2(i^)(l-/2(i^)))+ (34) i2f2Hil~f2ii^)) + h2-li-i{e-iy))fi{iy)il-f2iiy)) + 



((£-Z^)^ + (7l-72)')(/l(^)(l-/2(^)) + (l-/l(^))/2(i^))] 



Equation (^4|) gives the current noise at arbitrary tem- 
peratures and bias voltages. In the limit of small tem- 
peratures, it can be written as 



Si 



dvSi{v) , 



(35) 



Si{v) 



e^7i72 (e - + (7i - 72)^ 
7r2 [(e - z/)2 +72]2 



Note that the partial contribution Si{v)dv to the cur- 
rent noise (j3^) from an energy interval dv around some 
energy v is given by the standard general eqaation for the 
noise of a mesoscopic conductor - see, e.g., a, Si{v)dv — 
{e^dv / A-K^)D{1 — D), where the transparency D of the 
resonant-tunneling structure is energy-dependent: D = 
471 72 /((e — ly)^ + 7^)j and the width di/ of the energy 
interval acts as the bias energy eV. 

In the large-bias hmit, eq. (|5|) describes current 
noise associated with the inc pher mt sequential tunnel- 
ing through the resonant levelE3H'E2l: 



(72-7i+i(e-^))(l-/i('^))/2(i^))]. 
At T ^ 7, eV, this equation reduces to 



(37) 



527172 72 - 71 - *(£ - i^) 



[{e - J/)2 + 72]2 



Comparing eqs. (|27|), (p5D, and ( p8D for the low- 
temperature noise spectral densities, and the low- 
temperature limit of eq. (32) for the response coefficient, 
we can see that the partial contributions to these func- 
tions from tunneling of electrons with a certain energy v 
satisfy the generic relations of a linear quantum-limited 
detectorEJ: 



Si{y)SQ{v) = \Siq{v)\^ , IuySiq{v) 



-A(i.)/47r, 

(39) 



while, similarly to other types of detectors, ReS/Q char- 
acterizes the detector asymmetry. At small bias voltages, 
eV <^ 7, when the tunnel current is determined by tun- 
neling in the narrow energy range, the actual values of 
A and all spectral densities coincide with X{iy), Si{uX, 
etc. According to generic arguments presented in Ref.Eil, 
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eqs. ( |39D mean that the SET transistor in this regime 
is a quantum-hmited detector reaching the fundamental 
Hmit of the energy sensitivity e: 



27r 



[Si{iy)SQii^) - (Re5/Q(;.))2]i/2 ^ ri/2 . (40) 



Such a quantum detector also has an ultimate signal- 
to-noise ratio R for a continuous weak measurement of a 
qubit. If the qubit oscillates between its two basis states, 
the signal-to-noise ratio can be defined as the ratio of the 
amplitude of the oscillation line in the output spectrum 
of thcpdetector to the output noise. The ratio R can be 
shownEil to be closely related to the energy sensitivity 



R = 



A^^Si{u)Sq{u) 



For a symmetric detector with YLeSiQ{v) = (i.e., in 
the case of SET transistor, when 71 =72), the signal-to- 
noise ratio R can be written as (?i/e)^ and is limited by 
4.EJ This fundamental signal-to-noise ratio means that 
the detector does not add more dcphasing into the qubit 
dynamics than the minimum required by the|-acquisition 
of information about the state of the system£3'Ell 

The small bias condition necessary for the SET tran- 
sistor to operate as a quantum-limited detector with the 
energy sensitivity fi/2 limits the absolute value of the 
transistor output (the current /) . In order to see to what 
extent the transistor properties are degraded by increas- 
ing bias, we consider the case when the "collector" elec- 
trode absorbs all tunneling electrons, e — fi2 ^ J- De- 
tector characteristics of the transistor, most importantly 
energy sensitivity e, depend then on the position of the 
resonant energy level relative to the chemical potential of 
the "emitter" electrode characterized hy E = — e. In- 
tegrating eqs. (I?!), (11), (H), and (H) we find the noise 
spectral densities as functions of E: 



Si = ^??[(7? + iDil + tan-i z) - 27172 



l-l-z^ 



e^7i72f7r _i 
7r^7-' 2 



z = E/-f, (41) 



ReS'/g = (71 -72)5(5, -Im5_ 



IQ 



X 

47r 



1 



27r272 1 + z2 



Using eqs. (^ ) we obtain that the energy sensitivity 
e = 2n[SiSQ - (RbSiq)^]^/'^ / X of SET transistor in the 
resonant-tunneling regime is independent of the degree 
of the transistor asymmetry 71/72: 

^=5[(l + -^r(|+tan-iz)2-.T/^ (42) 



Equation (42) is plotted in Fig. |[ We see that the be- 
havior of e as a function of E changes from being nearly- 
constant for i? < to rapid increase with E at E > 0. 
This rapid increase reflects suppression of the response 
coefficient A ( |33| ) with increasing bias. Quantitatively, 
eq. (HI) gives for \E\ > 7: 



nE^/2-f^, E>0, 
l/VS, E<0. 



(43) 



This equation shows that in the regime of quantum tun- 
neling, E < 0, the energy sensitivity comes very close 
to the fundamental limit h/2. It does not quite reach it 
because of the energy dependence of the tunneling rate. 




FIG. 8. Low-temperature energy sensitivity ( ^^ of the 
SET transistor in the resonant-tunneling regime in units of 
h, a = e/Ti, as a function of the bias energy E characterizing 
position of the resonant level relative to the chemical potential 
of the emitter electrode. 

Figure ^ shows the signal-to-noise ratio R = 
A^/ (47r^ 5/ S'q) for continuous measurement of a qubit us- 
ing the SET transistor in the resonant-tunneling regime 
that is calculated from eqs. (^l|). In contrast to e, R de- 
pends on the transistor asymmetry, but the magnitude 
of variation of R with 71/72 is small, and regardless of 
71/72, the signal-to- noise ratio approaches 3 at large neg- 
ative E. In the symmetric case, 71/72 = 1, the signal- 
to-noise ratio is given by the general relation R — (?i/e)^ 
with energy sensitivity (^^. 

To conclude, we see that the SET transistor in the 
resonant-tunneling regime biased at the Coulomb block- 
ade threshold should operate at low temperatures as a 
quantum- limited detector. It has energy sensitivity Ti/VS 
and signal-to-noise ratio 3 for continuous weak measure- 
ment of a qubit. In this regime, transistor introduces the 
noise in the dynamics of the measured system that is only 
marginally larger than the fundamental noise required 
by the quantum mechanics of measurement. This means 
that for the optimally-biased transistor the time of the 
backaction dcphasing is equal to the measurement time. 
This conclusion differs from the results of the previows 
analysis of the SET transistor at large bias voltage£3. 
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where the transistor generates strong shot noise leading 
to the backaction dephasing that is much stronger than 
the fundamental minimum. 




FIG. 9. Signal-to-noise ratio R for continuous weak mea- 
surement of a qubit with the SET transistor in the reso- 
nant-tunneling regime. R is shown as a function of the bias 
energy E for various degrees of the transistor asymmetry. 

The author would like to acknowledge useful discus- 
sions with A.N. Korotkov. This work was supported by 
ARO and AFOSR. 



1 P.M.A. Dirac, Proc. R. Soc. A 114, 243 (1927). 

^ P.W. Anderson, in: "The many-body problem", ed. by E.R. 

Caianiello, (Academic Press, New York, 1964), vol. 2, p. 

113. 

^ "Quantum tunneling in condensed media", ed. by Yu. Ka- 
gan and A.J. Leggettt, (North-Holland, Amsterdam, 1992). 

* D.V. Averin and K.K. Likharev, in: Mesoscopic Phenom- 
ena m Solids, ed. by B.L. Altshuler, P.A. Lee, and R.A. 
Webb, (Elsevier, Amsterdam, 1991), p. 173. 

5 G. Schon and A.D. Zaikin, Phys. Rep. 198, 237 (1990). 
® "Single Charge Tunneling", ed. by H. Grabert and M. De- 
voret (Plenum, NY, 1992). 

D. Loss and D.P. DiVincenzo, Phys. Rev. A 57, 120 (1998). 

* Y. Nakamura, Yu.A. Pashkin, and J.S. Tsai, Nature 398, 
786 (1999). 

^ J.R. Fri edman, V. Patel, W . Chen, S.K. Tolpygo, and J.E. 

Lukens, |cond-mat/0004293| . 
^'^ J.E. Mooij, private communication. 

" D.V. Averin, A.B. Zorin, and K.K. Likharev, Zh. Eksp. 
Teor. Fiz. 88, 692 (1985) [Sov. Phys. JETP 61, 407]. 
D.T. Pegg and S.M. Barnett, J. Mod. Optics 44, 225 
(1997). 

" P.W. Anderson, J. Phys. Chem. Solids 11, 28 (1959). 
" K.A. Matveev and A.L Larkin, Phys. Rev. Lett. 78, 3749 
(1997). 



^® P.G. de Gennes, "Superconductivity of Metals and Alloys" 
(Addison- Wesley, 1989). 

W. Zerger, A.T. Dorsey, and M.P.A. Fisher, Phys. Rev. B, 

34, 6518 (1986). 
1^ S.M. Apenko, Phys. Lett. A 142, 277 (1989). 

D. Loss and K. Mullen, Phys. Rev. A 43, 2129 (1991). 
" Y.Nakamura, C.D. Chen, and J.S. Tsai, Phys. Rev. Lett. 

79, 2328 (1997). 

V. Bouchiat, D. Vion, P. Joyez, D. Esteve, and M.H. De- 
voret, J. of Supercond. 12, 789 (1999). 
■^^ D.J. Flees, S. Han, and J.E. Lukens, J. of Supercond. 12, 
813 (1999). 

D.V. Averin, J.R. Friedman, and J.E. Lukens, cond- 
mat/005081. 

D.V. Averin, J. Low Temp. Phys. 118, 781 (2000). 
Yu.N. Ovchinnikov and A. Schmidt, Phys. Rev. B 50, 6332 
(1994). 

■^^ R. Rouse, S. Han, and J.E. Lukens, Phys. Rev. Lett. 75, 
1614 (1995). 

^® P. Silvestrini, V.G. Palmieri, B. Ruggiero, and M. Russo, 
Phys. Rev. Lett. 79, 3046 (1997). 

K. Blum, "Density matrix theory and applications", 
(Plenum, NY, 1981). 

M. Biittiker, Phys. Rev. B 36, 3548 (1987). 
P. Lafarge, H. Pothier, E.R. Williams, D. Esteve, C. 
Urbina, and M.H. Devoret, Z. Phys. B 85, 327 (1991). 
^° P. Lafarge, P. Joyes, D. Esteve, C. Urbina, and M.H. De- 
voret, Phys. Rev. Lett 70, 994 (1993). 
D.V. Averin, Solid State Commun. 105, 659 (1998). 
Yu. Makhlin, G. Schon, and A. Shnirman, Nature 398, 305 

(1999) . 

^^L.B. loffe, V.B. Geshkenbein, M.V. Feigel'man, A.L. 
Fauchere, and G. Blatter, Nature 398, 679 (1999). 
J.E. Mooij, T.P. Orlando, L. Levitov, Lin Tian, Caspar H. 
van der Wal, and Seth Lloyd, Science 285, 1036 (1999). 
A. Blais and A. M. Zagoskin, Phys. Rev. A, 61 042308 

(2000) . 

C. Monroe, D.M. Meekhof, B.E. King, W.M. Itano, and 

D. J. Wineland, Phys. Rev. Lett. 75, 4714 (1995). 

^"^ I.L. Chuang, N.A. Gershenfeld, and M. Kubinec, Phys. 

Rev. Lett. 80, 3408 (1998). 
'■^^ R. Fazio, G.M. Palma, and J. Siewert, Phys. Rev. Lett. 83, 

5385 (1999). 

■^^ Yu. Makhlin, G. Schon, and A. Shnirman, J. Low Temp. 

Phys. 118, 751 (2000). 
*° D.V. Averin and K.K. Likharev, J. Low Temp. Phys. 62, 

345 (1986). 

*^ T.A. Fuhon and G.J. Dolan, Phys. Rev. Lett 59, 109 
(1987). 

*2 K.K. Likharev, IEEE Trans. Magn. 23, 1142 (1987). 

E. H. Visscher, S.M. Verbrugh, J. Lindeman, P. Hadley, and 
J.E. Mooij, AppL Phys. Lett. 66, 305 (1995). 

M.J. Yoo, T.A. Fulton, H.F. Hess, R.L. Willett, L.N. Dun- 

kleberger, R.J. Chichester, L.N. Pfeiffer, and K.W. West, 

Science 276, 579 (1997). 
*^ R.J. Schoelkopf, P. Wahlgren, A.A. Kozhevnikov, P. Dels- 

ing, and D.E. Prober, Science 280, 1238 (1998). 
"•^ B. Starmark, T. Henning, T. Claeson, P. Delsing, and A.N. 

Korotkov, J. Appl. Phys. 86, 2132 (1999). 
*^ V.A. Krupenin, D.E. Presnov, A.B. Zorin, and J. Niemeyer, 



12 



J. Low Temp. Phys. 118, 287 (2000). 

Yu.A. Pashkin, Y. Nakamura, and J.S. Tsai, Appl. Phys. 
Lett. 76, 2256 (2000). 

P.D. Dresselhaus, L. Ji, S. Han, J.E. Lukens, and K.K. 
Likharev, Phys. Rev. Lett 72, 3226 (1994). 
Y. Y. Wei, J. Weis, K. v. Khtzing, and K. Eberl, Phys. 
Rev. Lett 81, 1674 (1998). 

D.V. Averin, in: "Exploring the Quantum- Classical Fron- 
tier: Recent Advances m Macroscopic and Mesoscopic 
Quantum Pheno mena", Eds. J.R. F riedman and S. Han, 
to be published; cond-mat/0004364 , r-. 
D.V. Averin and Yu.V. Nazarov, in RefE, p. 217. 
'^'•^ A.N. Korotkov, D.V. Averin, K.K. Likharev, and S.A. 
Vasenko, in: "Single Electron Tunneling and Mesoscopic 
devices", Ed. by H. Koch and H. Liibbig, (Springer, Berlin, 
1992), p. 45. 

U. Hanke, Yu.M. Galperin, and K.A. Chao, Appl. Phys. 
Lett. 65, 1847 (1994); Erratum 65, 3416 (1994). 
Yu.V. Nazarov, J. Low Temp. Phys. 90, 77 (1993). 
P. Lafarge and D. Esteve, Phys. Rev. B 48, 14309 (1993). 
" D.V. Averin, Physica B, 194-196, 979 (1994). 

H. Schoeller and G. Schon, Phys. Rev. B 50, 18436 (1994). 
L.P. Kouwenhoven, CM. Marcus, P.L. McEuen, S. 
Tarucha, R.M. Westervelt, and N.S. Wingreen, in: "Meso- 
scopic Electron Transport", Ed. by L.L. Sohn, L.P. 
Kouwenhoven, and G. Schon, (Kluwer, 1997), p. 105. 
D.V. Averin, J. Appl. Phys. 73, 2593 (1£^)- 
The charge sensitivity discussed in Ref.tJ was expressed 
through the energy sensitivity 5e of the tunnel current to 
variations of the position of the resonant level. Note that 
this energy sensitivity is different from the more standard 
energy sensitivity e (WG) studied in this work. 
L.Y. Chen and C.S. Ting, Phys. Rev. B 43, 4534 (1991). 
J.H. Davies, P. Hyldgaard, S. Hershfield, and J.W. Wilkins, 
Phys. Rev. B 46, 9620 (1992). 



Ya.M. Blanter and M. Buttiker, |cond-mat/991015^. 
A.N. Korotkov and D.V. Averin, ^ond-mat 70002203 . 
'''^ A.N. Korotkov, Phys. Rev. B 60, 5737 (1999). 



67 



A. Shnirman and G. Schon, Phys. Rev. B 57, 15400 (1998). 



13 



